An Immune Gene Expression Risk Score for Distant Metastases after Radiotherapy for Cervical Cancer.

PURPOSE: To develop an immune-based gene expression risk score to identify patients with cervical cancer at increased risk of distant metastases (DM). EXPERIMENTAL DESIGN: Tumor biopsies were obtained from 81 patients prior to chemoradiotherapy. Whole-transcriptome RNA sequencing was performed (Illumina NextSeq500). Beginning with 4,723 immune-related genes, a 55-gene risk score for DM was derived using Cox modeling and principal component analysis. It was validated in independent cohorts of 274 patients treated at the Norwegian Radium Hospital (NRH) and 206 patients from The Cancer Genome Atlas (TCGA). RESULTS: The risk score was predictive of DM (HR, 2.7; P < 0.0001) and lower cause-specific survival (CSS) by univariate analysis (HR, 2.0; P = 0.0003) and multivariate analysis adjusted for clinical factors (DM HR, 3.0; P < 0.0001; CSS HR, 2.2; P = 0.0004). The risk score predicted DM (HR, 1.4; P = 0.05) and CSS (HR, 1.48; P = 0.013) in the NRH cohort and CSS (HR, 1.4; P = 0.03) in TCGA cohort. Higher risk scores were associated with lower CIBERSORT estimates of tumor-infiltrating immune cells, including CD8 T cells and M1 and M2 macrophages (all P < 0.001). Higher risk scores were associated with lower expression (all P < 0.001) of important chemokines (CXCL12, CXCR4), IFN-regulated genes (IRF1, STAT1, IDO1), and immune checkpoint regulators (PD-1, PD-L1, CTLA-4). CONCLUSIONS: The immune metastatic risk score addresses important challenges in the treatment of cervical cancer-identifying patients at high risk of DM after radiotherapy. The findings of this study indicate that high tumor mutational burden and a "cold," immune-excluded tumor microenvironment influence distant metastatic recurrence. Further validation of the risk score is needed.

ALDH1A1 drives prostate cancer metastases and radioresistance by interplay with AR- and RAR-dependent transcription.

Rationale: Current therapies for metastatic osseous disease frequently fail to provide a durable treatment response. To date, there are only limited therapeutic options for metastatic prostate cancer, the mechanisms that drive the survival of metastasis-initiating cells are poorly characterized, and reliable prognostic markers are missing. A high aldehyde dehydrogenase (ALDH) activity has been long considered a marker of cancer stem cells (CSC). Our study characterized a differential role of ALDH1A1 and ALDH1A3 genes as regulators of prostate cancer progression and metastatic growth. Methods: By genetic silencing of ALDH1A1 and ALDH1A3 in vitro, in xenografted zebrafish and murine models, and by comparative immunohistochemical analyses of benign, primary tumor, and metastatic specimens from patients with prostate cancer, we demonstrated that ALDH1A1 and ALDH1A3 maintain the CSC phenotype and radioresistance and regulate bone metastasis-initiating cells. We have validated ALDH1A1 and ALDH1A3 as potential biomarkers of clinical outcomes in the independent cohorts of patients with PCa. Furthermore, by RNAseq, chromatin immunoprecipitation (ChIP), and biostatistics analyses, we suggested the molecular mechanisms explaining the role of ALDH1A1 in PCa progression. Results: We found that aldehyde dehydrogenase protein ALDH1A1 positively regulates tumor cell survival in circulation, extravasation, and metastatic dissemination, whereas ALDH1A3 plays the opposite role. ALDH1A1 and ALDH1A3 are differentially expressed in metastatic tumors of patients with prostate cancer, and their expression levels oppositely correlate with clinical outcomes. Prostate cancer progression is associated with the increasing interplay of ALDH1A1 with androgen receptor (AR) and retinoid receptor (RAR) transcriptional programs. Polo-like kinase 3 (PLK3) was identified as a transcriptional target oppositely regulated by ALDH1A1 and ALDH1A3 genes in RAR and AR-dependent manner. PLK3 contributes to the control of prostate cancer cell proliferation, migration, DNA repair, and radioresistance. ALDH1A1 gain in prostate cancer bone metastases is associated with high PLK3 expression. Conclusion: This report provides the first evidence that ALDH1A1 and PLK3 could serve as biomarkers to predict metastatic dissemination and radiotherapy resistance in patients with prostate cancer and could be potential therapeutic targets to eliminate metastasis-initiating and radioresistant tumor cell populations.

Integrated analysis of cervical squamous cell carcinoma cohorts from three continents reveals conserved subtypes of prognostic significance.

Human papillomavirus (HPV)-associated cervical cancer is a leading cause of cancer deaths in women. Here we present an integrated multi-omic analysis of 643 cervical squamous cell carcinomas (CSCC, the most common histological variant of cervical cancer), representing patient populations from the USA, Europe and Sub-Saharan Africa and identify two CSCC subtypes (C1 and C2) with differing prognosis. C1 and C2 tumours can be driven by either of the two most common HPV types in cervical cancer (16 and 18) and while HPV16 and HPV18 are overrepresented among C1 and C2 tumours respectively, the prognostic difference between groups is not due to HPV type. C2 tumours, which comprise approximately 20% of CSCCs across these cohorts, display distinct genomic alterations, including loss or mutation of the STK11 tumour suppressor gene, increased expression of several immune checkpoint genes and differences in the tumour immune microenvironment that may explain the shorter survival associated with this group. In conclusion, we identify two therapy-relevant CSCC subtypes that share the same defining characteristics across three geographically diverse cohorts.

miR-200a/b/-429 downregulation is a candidate biomarker of tumor radioresistance and independent of hypoxia in locally advanced cervical cancer.

Many patients with locally advanced cervical cancer experience recurrence within the radiation field after chemoradiotherapy. Biomarkers of tumor radioresistance are required to identify patients in need of intensified treatment. Here, the biomarker potential of miR-200 family members was investigated in this disease. Also, involvement of tumor hypoxia in the radioresistance mechanism was determined, using a previously defined 6-gene hypoxia classifier. miR-200 expression was measured in pretreatment tumor biopsies of an explorative cohort (n = 90) and validation cohort 1 (n = 110) by RNA sequencing. Publicly available miR-200 data of 79 patients were included for the validation of prognostic significance. A score based on expression of the miR-200a/b/-429 (miR-200a, miR-200b, and miR-429) cluster showed prognostic significance in all cohorts. The score was significant in multivariate analysis of central pelvic recurrence. No association with distant recurrence or hypoxia status was found. Potential miRNA target genes were identified from gene expression profiles and showed enrichment of genes in extracellular matrix organization and cell adhesion. miR-200a/b/-429 overexpression had a pronounced radiosensitizing effect in tumor xenografts, whereas the effect was minor in vitro. In conclusion, miR-200a/b/-429 downregulation is a candidate biomarker of central pelvic recurrence and seems to predict cell adhesion-mediated tumor radioresistance independent of clinical markers and hypoxia.

MRI Distinguishes Tumor Hypoxia Levels of Different Prognostic and Biological Significance in Cervical Cancer.

Tumor hypoxia levels range from mild to severe and have different biological and therapeutical consequences but are not easily assessable in patients. Here we present a method based on diagnostic dynamic contrast enhanced (DCE) MRI that reflects a continuous range of hypoxia levels in patients with tumors of cervical cancer. Hypoxia images were generated using an established approach based on pixel-wise combination of DCE-MRI parameters ν e and K trans, representing oxygen consumption and supply, respectively. Using two tumor models, an algorithm to retrieve surrogate measures of hypoxia levels from the images was developed and validated by comparing the MRI-defined levels with hypoxia levels reflected in pimonidazole-stained histologic sections. An additional indicator of hypoxia levels in patient tumors was established on the basis of expression of nine hypoxia-responsive genes; a strong correlation was found between these indicator values and MRI-defined hypoxia levels in 63 patients. Chemoradiotherapy outcome of 74 patients was most strongly predicted by moderate hypoxia levels, whereas more severe or milder levels were less predictive. By combining gene expression profiles and MRI-defined hypoxia levels in cancer hallmark analysis, we identified a distribution of levels associated with each hallmark; oxidative phosphorylation and G2-M checkpoint were associated with moderate hypoxia, epithelial-to-mesenchymal transition, and inflammatory responses with significantly more severe levels. At the mildest levels, IFN response hallmarks together with HIF1A protein expression by IHC appeared significant. Thus, our method visualizes the distribution of hypoxia levels within patient tumors and has potential to distinguish levels of different prognostic and biological significance. SIGNIFICANCE: These findings present an approach to image a continuous range of hypoxia levels in tumors and demonstrate the combination of imaging with molecular data to better understand the biology behind these different levels.

Combining imaging- and gene-based hypoxia biomarkers in cervical cancer improves prediction of chemoradiotherapy failure independent of intratumour heterogeneity.

BACKGROUND: Emerging biomarkers from medical imaging or molecular characterization of tumour biopsies open up for combining the two and exploiting their synergy in treatment planning of cancer patients. We generated a paired data set of imaging- and gene-based hypoxia biomarkers in cervical cancer, appraised the influence of intratumour heterogeneity in patient classification, and investigated the benefit of combining the methodologies in prediction of chemoradiotherapy failure. METHODS: Hypoxic fraction from dynamic contrast enhanced (DCE)-MR images and an expression signature of six hypoxia-responsive genes were assessed as imaging- and gene-based biomarker, respectively in 118 patients. FINDINGS: Dichotomous biomarker cutoff to yield similar hypoxia status by imaging and genes was defined in 41 patients, and the association was validated in the remaining 77 patients. The two biomarkers classified 75% of 118 patients with the same hypoxia status, and inconsistent classification was not related to imaging-defined intratumour heterogeneity in hypoxia. Gene-based hypoxia was independent on tumour cell fraction in the biopsies and showed minor heterogeneity across multiple samples in 9 tumours. Combining imaging- and gene-based classification gave a significantly better prediction of PFS than one biomarker alone. A combined dichotomous biomarker optimized in 77 patients showed a large separation in PFS between more and less hypoxic tumours, and separated the remaining 41 patients with different PFS. The combined biomarker showed prognostic value together with tumour stage in multivariate analysis. INTERPRETATION: Combining imaging- and gene-based biomarkers may enable more precise and informative assessment of hypoxia-related chemoradiotherapy resistance in cervical cancer. FUNDING: Norwegian Cancer Society, South-Eastern Norway Regional Health Authority, and Norwegian Research Council.

Identification and Validation of Reference Genes for RT-qPCR Studies of Hypoxia in Squamous Cervical Cancer Patients.

Hypoxia is an adverse factor in cervical cancer, and hypoxia-related gene expression could be a powerful biomarker for identifying the aggressive hypoxic tumors. Reverse transcription quantitative PCR (RT-qPCR) is a valuable method for gene expression studies, but suitable reference genes for data normalization that are independent of hypoxia status and clinical parameters of cervical tumors are lacking. In the present work, we aimed to identify reference genes for RT-qPCR studies of hypoxia in squamous cervical cancer. From 422 candidate reference genes selected from the literature, we used Illumina array-based expression profiles to identify 182 genes not affected by hypoxia in cervical cancer, i.e. genes regulated by hypoxia in eight cervical cancer cell lines or correlating with the hypoxia-associated dynamic contrast-enhanced magnetic resonance imaging parameter ABrix in 42 patients, were excluded. Among the 182 genes, nine candidates (CHCHD1, GNB2L1, IPO8, LASP1, RPL27A, RPS12, SOD1, SRSF9, TMBIM6) that were not associated with tumor volume, stage, lymph node involvement or disease progression in array data of 150 patients, were selected for further testing by RT-qPCR. geNorm and NormFinder analyses of RT-qPCR data of 74 patients identified CHCHD1, SRSF9 and TMBIM6 as the optimal set of reference genes, with stable expression both overall and across patient subgroups with different hypoxia status (ABrix) and clinical parameters. The suitability of the three reference genes were validated in studies of the hypoxia-induced genes DDIT3, ERO1A, and STC2. After normalization, the RT-qPCR data of these genes showed a significant correlation with Illumina expression (P<0.001, n = 74) and ABrix (P<0.05, n = 32), and the STC2 data were associated with clinical outcome, in accordance with the Illumina data. Thus, CHCHD1, SRSF9 and TMBIM6 seem to be suitable reference genes for studying hypoxia-related gene expression in squamous cervical cancer samples by RT-qPCR. Moreover, STC2 is a promising prognostic hypoxia biomarker in cervical cancer.

Integrative Analysis of DCE-MRI and Gene Expression Profiles in Construction of a Gene Classifier for Assessment of Hypoxia-Related Risk of Chemoradiotherapy Failure in Cervical Cancer.

PURPOSE: A 31-gene expression signature reflected in dynamic contrast enhanced (DCE)-MR images and correlated with hypoxia-related aggressiveness in cervical cancer was identified in previous work. We here aimed to construct a dichotomous classifier with key signature genes and a predefined classification threshold that separated cervical cancer patients into a more and less hypoxic group with different outcome to chemoradiotherapy. EXPERIMENTAL DESIGN: A training cohort of 42 patients and two independent cohorts of 108 and 131 patients were included. Gene expression data were generated from tumor biopsies by two Illumina array generations (WG-6, HT-12). Technical transfer of the classifier to a reverse transcription quantitative PCR (RT-qPCR) platform was performed for 74 patients. The amplitude ABrix in the Brix pharmacokinetic model was extracted from DCE-MR images of 64 patients and used as an indicator of hypoxia. RESULTS: Classifier candidates were constructed by integrative analysis of ABrix and gene expression profiles in the training cohort and evaluated by a leave-one-out cross-validation approach. On the basis of their ability to separate patients correctly according to hypoxia status, a 6-gene classifier was identified. The classifier separated the patients into two groups with different progression-free survival probability. The robustness of the classifier was demonstrated by successful validation of hypoxia association and prognostic value across cohorts, array generations, and assay platforms. The prognostic value was independent of existing clinical markers, regardless of clinical endpoints. CONCLUSIONS: A robust DCE-MRI-associated gene classifier has been constructed that may be used to achieve an early indication of patients' risk of hypoxia-related chemoradiotherapy failure. Clin Cancer Res; 22(16); 4067-76. ©2016 AACR.

Interplay between promoter methylation and chromosomal loss in gene silencing at 3p11-p14 in cervical cancer.

Loss of 3p11-p14 is a frequent event in epithelial cancer and a candidate prognostic biomarker in cervical cancer. In addition to loss, promoter methylation can participate in gene silencing and promote tumor aggressiveness. We have performed a complete mapping of promoter methylation at 3p11-p14 in two independent cohorts of cervical cancer patients (n = 149, n = 121), using Illumina 450K methylation arrays. The aim was to investigate whether hyperm-ethylation was frequent and could contribute to gene silencing and disease aggressiveness either alone or combined with loss. By comparing the methylation level of individual CpG sites with corresponding data of normal cervical tissue, 26 out of 41 genes were found to be hypermethylated in both cohorts. The frequency of patients with hypermethylation of these genes was found to be higher at tumor stages of 3 and 4 than in stage 1 tumors. Seventeen of the 26 genes were transcriptionally downregulated in cancer compared to normal tissue, whereof 6 genes showed a significant correlation between methylation and expression. Integrated analysis of methylation, gene dosage, and expression of the 26 hypermethylated genes identified 3 regulation patterns encompassing 8 hypermethylated genes; a methylation driven pattern (C3orf14, GPR27, ZNF717), a gene dosage driven pattern (THOC7, PSMD6), and a combined methylation and gene dosage driven pattern (FHIT, ADAMTS9, LRIG1). In survival analysis, patients with both hypermethylation and loss of LRIG1 had a worse outcome compared to those harboring only hypermethylation or none of the events. C3orf14 emerged as a novel methylation regulated suppressor gene, for which knockdown was found to promote invasive growth in human papilloma virus (HPV)-transformed keratinocytes. In conclusion, hypermethylation at 3p11-p14 is common in cervical cancer and may exert a selection pressure during carcinogenesis alone or combined with loss. Information on both events could lead to improved prognostic markers.

The duration of gastrin treatment affects global gene expression and molecular responses involved in ER stress and anti-apoptosis.

BACKGROUND: How cells decipher the duration of an external signal into different transcriptional outcomes is poorly understood. The hormone gastrin can promote a variety of cellular responses including proliferation, differentiation, migration and anti-apoptosis. While gastrin in normal concentrations has important physiological functions in the gastrointestine, prolonged high levels of gastrin (hypergastrinemia) is related to pathophysiological processes. RESULTS: We have used genome-wide microarray time series analysis and molecular studies to identify genes that are affected by the duration of gastrin treatment in adenocarcinoma cells. Among 403 genes differentially regulated in transiently (gastrin removed after 1 h) versus sustained (gastrin present for 14 h) treated cells, 259 genes upregulated by sustained gastrin treatment compared to untreated controls were expressed at lower levels in the transient mode. The difference was subtle for early genes like Junb and c-Fos, but substantial for delayed and late genes. Inhibition of protein synthesis by cycloheximide was used to distinguish between primary and secondary gastrin regulated genes. The majority of gastrin upregulated genes lower expressed in transiently treated cells were primary genes induced independently of de novo protein synthesis. This indicates that the duration effect of gastrin treatment is mainly mediated via post-translational signalling events, while a smaller fraction of the differentially expressed genes are regulated downstream of primary transcriptional events. Indeed, sustained gastrin treatment specifically induced prolonged ERK1/2 activation and elevated levels of the AP-1 subunit protein JUNB. Enrichment analyses of the differentially expressed genes suggested that endoplasmic reticulum (ER) stress and survival is affected by the duration of gastrin treatment. Sustained treatment exerted an anti-apoptotic effect on serum starvation-induced apoptosis via a PKC-dependent mechanism. In accordance with this, only sustained treatment induced anti-apoptotic genes like Clu, Selm and Mcl1, while the pro-apoptotic gene Casp2 was more highly expressed in transiently treated cells. Knockdown studies showed that JUNB is involved in sustained gastrin induced expression of the UPR/ER stress related genes Atf4, Herpud1 and Chac1. CONCLUSION: The duration of gastrin treatment affects both intracellular signalling mechanisms and gene expression, and ERK1/2 and AP-1 seem to play a role in converting different durations of gastrin treatment into distinct cellular responses.

Gastrin-induced proliferation involves MEK partner 1 (MP1).

The peptide hormone gastrin is an important factor for the maintenance and homeostasis of the gastric mucosa. We show that gastrin stimulates proliferation in a dose-dependent manner in the human gastric adenocarcinoma cell line AGS-GR. Furthermore, we demonstrate that the MAPK scaffold protein MEK partner 1 (MP1) is important for gastrin-induced phosphorylation of ERK1 and ERK2 and that MP1 promotes gastrin-induced proliferation of AGS-GR cells. Our results suggest a role of MP1 in gastrin-induced cellular responses involved in proliferation and homeostasis of the gastric mucosa.